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[57] ABSTRACT 

A computer implemented risk-management system sched- 
ules the generating units of an electric utility while taking 
into consideration power trading with other utilities and the 
stochastic load on the utility system. The system provides 
the user with a tool that generates multiple load forecasts and 
allows the user to vary the fuel price between the different 
scenarios and the different periods of the planning horizon. 
The tool allows the user to model accurately the uncertain 
trading transactions and the changing fuel prices to meet the 
electric demand of customers al a minimal cost while 
making the maximum profit possible from power trading. 
The tool also allows the user to apply any set of linear 
constraints to fuels. A mathematical model of the problem is 
solved to provide the status of each generator at each time 
period of the planning horizon under each given scenario, 
the load on each generator during each period in which it is 
operating, an optimal fuel mix for each generating unit, and 
the prices for purchasing and selling power in the periods of 
the planning horizon. 

7 Claims 13 Drawing Sheets. 
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FORECASTED TRADES FOR OCTOBER 15, 1996 


TYPE 


MWH 


$/MWH 


TIME 


PROB 


BUY 


50 


15.00 


OFF-PEAK 


25% 


BUY 


100 


16.00 


OFF-PEAK 


25% 


BUY 


150 


17.00 


ON-PEAK 


30% 


BUY 


175 


19.00 


ON-PEAK 


40% 


SELL 


100 


22.00 


OFF-PEAK 


20% 







. UNIT 


TYPE 


MIN RUN 


MIN DOWN 


MIN LOAD 


MAX LOAD 


FUEL 


$/MMBTUI 


!! ANN 1 


R 
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1 


200 


500 


GAS 1 


2.00 


. ANN 2 


R 


1 


1 


200 


500 


GAS 2 


2.10 


. FRAN 1 
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4 


8 


50 


110 


COAL 1 


1.3 


. FRAN 2 
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8 


50 


110 


COAL 2 


1.4 


. SAM 


C 
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4 


25 


80 


COAL 3 


1.6 













FORECASTED LOAD 
FOR OCTOBER 15, 
1996 


HOUR 


LOAD 


00;00 


10495 


01:00 


11135 


02:00 


11653 


03:00 


11829 


04:00 


11995 



2338 



2490 



2787 



V 



RISK MANAGEMENT SYSTEM 



SCHEDULE 
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01:00 


ANN 1 
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270 


ANN 2 


231 


270 


FRAN 1 
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50 


FRAN 2 


50 


50 
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FUEL 


00:00 


01:00 


02:00 


03:00 


04:00 


GAS 1 


12 
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GAS 2 
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COAL 1 


23 


24 


24 


26 


28 


COAL 2 
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10 


12 


14 


15 


COAL 3 
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0 
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INPUT LOAD FORECASTS AND GENERATOR 
PROPERTIES (INCLUDING COST CURVES). 



USE LAGRANGIAN RELAXATION TO DECOMPOSE 
THE PROBLEM INTO SINGLE-GENERATOR 
SUBPROBLEMS. 



INITIALIZE THE LAGRANGIAN MULTIPLIERS. 



SOLVE THE (STOCHASTIC) DYNAMIC PROGRAM 
OF EACH GENERATOR. 




UPDATE THE LAGRANGE MULTIPLIERS. 



OUTPUT GENERATION LEVELS UNDER 
DIFFERENT SCENARIOS. 




41 



42 



43 



44 



FIG.4 

PRIOR ART 
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INPUT LOAD FORECASTS AND GENERATOR 
PROPE RTIES (INCLU DING COST CURVES). 



USE LAGRANGIAN RELAXATION (PROGRESSIVE 
HEDGING) TO RELAX THE BUNDLE CONSTRAINTS, 
EACH SCENARIO HAS ITS OWN PROBLEM. 



I 



INITIALIZE THE PROGRESSIVE HEDGING 
PENALTIES. 



50 



51 



52 



USE LAGRANGIAN RELAXATION TO DECOMPOSE 
THE PROBLEM INTO SINGLE-GENERATOR 
SUBPROBLEMS. 



INITIALIZE THE LAGRANGIAN MULTIPLIERS 



SOLVE THE (STOCHASTIC) DYNAMIC PROGRAM 
OF EACH GENERATOR. 




UPDATE PROGRESSIV E HEDGING PENALTIES. 

1 



£ 



OUTPUT GENERATION LEVELS UNDER 
DIFFERENT SCENARIOS. 



FIG.5 
PRIOR ART 




TRADITIONAL 
UNIT 

COMMITMENT 
MODEL 
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INPUT LOAD FORECAST, TRADING FORECAST, 
AND GENERATOR AND FUEL PROPERTIES. 



71 



CALL TREE 
I 



]_/~72 



CALL BUNDLE IF NEEDED 





APPROXIMATE THE COST FUNCTION OF EACH 
GENERATOR BY USING THE CHEAPEST FUEL 










CALL SUCP TO OBTAIN LOAD ON GENERATORS 
UNDER ANY SCENARIO. 



p-73 
74 



75 



GIVEN THE LOAD ON EACH GENERATOR, CALL 
FUEL TO DETERMINE AN OPTIMAL FUEL 
ALLOCATION. 




COMPUTE A NEW COST FUNCTION FOR EACH 
GENERATOR. 



76 



OUTPUT GENERATION LEVELS, PRICING 
INFORMATION. AND FUEL USAGE. 




FIG.7 
PRIOR ART 
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FIG. 10 



SET THE PREDECESSOR OF SCENARIO 1 TO 1 
AND ITS SPUT TIME TO 1. LET s=2. 



COMPARE THE LOAD OF s WITH THE LOAD OF 
1 TO OBTAIN THE SPLIT TIME OF s. SET 
THE PREDECESSOR OF s TO 1. SET s=2. 




COMPARE THE LOADS OF s AND s TO OBTAIN 
A NEW SPLIT TIME FOR s. LET s BE THE 
THE PREDECESSOR OF s. 



106 



INCREMENT s* BY 1 



INCREMENT s BY 1 



NO 




101 



102 



107 
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12 1 ^ \ INITIALIZE THE LAGRANGE MULTIPLIERS 

i 

122 USE THE BUNDLE-TRUST REGION METHOD TO 
SOLVE THE DUAL PROBLEM. 



i 



12 3 DETERMINE A FEASIBLE SOLUTION FOR THE 
PRIMAL PROBLEM. 



i 



12 4 ^ SOLVE AN ECONOMIC POWER DISPATCH 




FIG. 12 
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SET beto(s) TO T FOR ALL SCENARIOS s. 
SET s=S. 



1401 



SOLVE THE DYNAMIC PROGRAM OF SCENARIO 
s BY MOVING BACKWARD FROM beta(s) TO 
tau(s). 



1402 



SOLVE THE DYNAMIC PROGRAM OF sigmo(s) 
BY MOVING BACKWARD FROM beta(s) TO 
tau(s). INCREASE p(sigma(s)) BY p(s). 






SET beta(sigma(s)) TO tau(s)-1. 






DECREASE s BY 1. 




' /-1406 



NO 




GIVEN THE INITIAL STATUS OF THE GENERATOR, 
MOVE FORWARD TO TRACE THE OPTIMAL 
DECISIONS. 



1403 



1404 



1405 



1407 



(end) 



FIG. 14 
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SET k=1, s=2, t=1, i=1, AND ]=1. 



1501 



SET B(k, l(y(s, i, t, j))} = 1 AND 
B(k, l(sigma(s), i, t, jj) ~ 1. 
INCREMENT k BY 1. 



1502 



NO 



NO 



NO 



NO 



FIG. 15 



INCREMENT j BY 1 



INCREMENT t BY j 



(end) 
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RISK MANAGEMENT SYSTEM FOR per MWH or the "unit heat-rate curve". Note that the unit 

ELECTRIC UTILITIES heat-rate curve shown in FIG. IB starts at a maximum rate 

(BTU/MWH) and declines as the load increases. This 

BACKGROUND OF THE INVENTION reflects the fact that the unit's efficiency increases as we 

1. Field of the Invention 5 increase the load towards the middle of the operating range. 
Hie present invention generally relates to computer sys- n * efficiency iisuallydecUnes after a certain load, x To 

terns for scheduling the generating units of electric utilities * ranslate th u e curv f of nGS * 1A and * B mt0 a relationship 

and, more particularly, to a computer implemented process load 0D a 1??^' and me u c °5 one can r 

which, given (1) a set of load forecast on the electric SS??a H°?S ^ } ^ 

system, (2) a prediction of trading transactions that may take 10 FIGS * 1Aand 1B - 

place within the next few days, (3) a description of the The other classes of generators, nuclear and hydroelectric, 

physical properties of the generating units, (4) an estimate of use similar concepts to generate electricity. In the case of 

spot prices, and (5) fuel prices associated with load nuclear units, a radioactive material is used to run the 

forecasts, provides generation levels, fuel usage and prices reactor; while in the case of hydroelectric units, a water 

at which power can be traded. 15 stream is used to push the turbines. Note that thermal and 

2. Background Description nuclear units require a lead-time (two hours or more) to start 

ci . • • i -* c i* ir t i generating power. Given that there might be an urgent need 

Electricity is an essential part of our lives. If we took a f **'„.. . . ..^ & & * 

_ * * #u* i r .u ■ *u * u . • * r for power, utilities maintain quick-start generators for 

moment to think of the services that would not exist if i • * -m. \ . i . 

, ^- . , , , . j A1 unusual circumstances. These units use natural gas to heat 

electricity were not here, we would be amazed. Almost - n . 4 ,. , . . , , . . 

every aspect of our modern lives involves electric power; 20 n °* ™ ter - wh , lc * K then 10 P uah . 1 smM P"*? 6 - 

from light bulbs and television sets to hospitals and auto pearly these units have a smaller generating capacity than 

■a * : Am, u i* u • u thermal units and are not as efficient. However, they have a 

industries. Although we are used to having power whenever , . . ^ ✓ • v , 

we need it, the processes and systems involved in delivering ver y short star V u P ' lm t e ( ten mmu ,f or k less ) because the 

electricity require careful planning and sophisticated math „ energy required to heat air is less than that needed to heat 

ematical and forecasting models. 25 water. The quick-start generators are similar m their opera- 

° tionc in lot pnntnuc 



Lately, due to rising costs of energy and discrepancies in 



tions to jet engines. 



its price among different regions of the country, the legal t In terms °. f °Pf rall0ns > generators can be classified into 

framework surrounding the electric-power industry has three categories: must-run units , cyclers .and peakers . 

changed. This change has opened the door for utilities to 30 M^'"™ 1 ™ lts need . to ^ ™ am & aU ^ e . nme - ™ ere | s a 

compete with each other and against independent suppliers ™mber of reasons for this. One example is that of nuclear 

regardless of their geographic location. Although this units The time and energy required to start them are costly 

change will benefit the consumer, utilities are going to face which makes switching them on and off economically 

a highly unpredictable market and will need to make tough feasible Another example is when a utility goes into a 

decisions regarding power generation and delivery. 35 contractual agreement with a generating plant which results 

~ , c m continual running for a certain penod. 

Depending on the fuel used, electric-power generators of or 

a utility may be categorized into thermal, nuclear, and ^ second category of generators is cyclers. These are 

hydroelectric. In general, a thermal unit uses a boiler that units that are switched on and off depending on need. To 

generates steam to drive a turbine. The turbine motion is avoid damaging these units, once a unit has been switched 

then utilized to generate electric power. After cooling down, 40 on » il must remain running for a minimal time duration (four 

condensed steam (water) is pumped back into the boiler so to fort y-eight hours). A similar constraint holds when a 

that it can be heated again. Part of the power generated is generator is switched off; it has to be given enough time to 

used to drive a cooling fan and to run the system's pump. cog1 down ( four t0 forty-eight hours). In general, starting up 

Both coal and oil can be used to fuel the boiler. To start a a c y cler ma y take two to twenty-four hours depending on the 

thermal generator quickly, natural gas is often used; then the 45 unit * If P ower » needed ur g entlv ' a utim y uses P eakill g unit ** 

system switches to the standard fuel. Note that the opera- are the quick-start units mentioned above. They are 

tions of a thermal unit are similar to those of a locomotive. cosll y but necessary to maintain a reliable system that can 

The output of a thermal unit is a function of the heat deliver P ower under an y circumstances, 

generated in the boiler. Given that there is a limit to the The generated electric -power is carried across transmis- 

amount of fuel that can be burnt during one hour, a thermal 50 sion lines to the different points of usage or delivery points, 

generator has a maximum generation capacity, G. Generat- Examples of delivery points are a city or a large factory, 

ing enough steam to drive the turbine requires a minimal Whenever electric power is transmitted, part of it is lost. To 

amount of fuel which depends on the generating-unit's maintain a reliable service, a utility must generate, at any 

physical properties. Hence, we expect the minimum output moment, enough electricity to meet the demand of its 

of a generator, g, to be significantly greater than zero. To 55 customers plus losses incurred in power transmission. Note 

start a generator, the boiler must first be heated to the level that exceeding customers' need increases the cost of oper- 

at which it starts producing steam. The cost of the fuel used ating the system and may cause power outages due to 

in heating the boiler is called the "start-up" cost. Note that overloading some power lines. Given that power usage at 

this cost is not associated with any electric-power output. any delivery point fluctuates continually, the problem of 

FIGS. 1A and IB present the relationship between the 60 controlling the load on each generator while meeting the 

heat input to a generator and its electric-power output. The demand requires sophisticated tools, 

heat input is measured in British Thermal Units or BTU. The To increase the reliability of an electric-power system, 

curve in FIG. 1A representing the electric-power output transmission lines of different utilities are connected with 

(MWH) is usually a second-degree polynomial in heat input. each other at certain points. These connections allow one 

As expected, this function is increasing; i.e., as we burn 65 utility to borrow power from another utility if needed. For 

more fuel we expect to generate more electricity. One can example, a generator may fail suddenly resulting in shortage 

use the curve in FIG. 1A to compute the average heat input of generation. Being connected to neighboring utilities, the 
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burden of generation loss is shared with others. Of course, pricing scheme. Deregulation is changing this picture by 
the utility that suffered a generator loss must start another un-bundling the electric power into generation and trans- 
unit as quickly as possible or take action to compensate for mission. One will pay a production cost and a transmission 
this loss. fee. There will be several power suppliers from whom 
To meet customers' demand, a planner predicts the elec- 5 electric power may be purchased. Suppliers may have dif- 
tric load on the system for the following week. The forecast ferent pricing mechanisms. For example, there might be a 
is given in the form of total load per hour. That is, it gives discount for using power off-peak periods or for signing a 
the expected load at each hour for the next 168 hours. long-term contract with the supplier. Power producers will 
Electric utilities apply advanced statistical models to the compete with each other to minimize their costs so that they 
weather forecast for each region they serve, load data from ™ can sell their product to more customers and maximize their 
previous years, and load data from last week; and use profit. 

intuition to produce an accurate forecast. Using the load On the other hand, power transmission will remain regu- 

forecast, a utility schedules its generating units for the next lated for the time being. The reason for that is to maintain a 

168 hours. By scheduling, what is meant is deciding when reliable system. The transmission lines in each state or 

each generator should be on or off during the week. The 15 region will be controlled by an independent entity called 

schedule also indicates the load on each generator in the Independent System Operator or ISO. One of the ISO 

system at each hour. To minimize the cost of operating the responsibilities is to settle financially with the parties 

system, a utility attempts to find a schedule that meets the involved in transmitting electric power. The transmission 

load at a minimal cost. It is important to note that the cost depends on the proximity of the supplier and the 

resulting schedule must take into consideration all con- 20 congestion of the transmission lines as well as other opera- 

straints imposed by the system, such as minimum up-time tional factors. To maintain a reliable system, ISO announces, 

and minimum down-time, which makes the scheduling task 24 hours in advance, the load forecast on the system and asks 

notoriously hard. The problem of rinding a minimal cost interested suppliers to submit bids. The ISO then holds an 

schedule for the generating units is known as the "unit auction to determine which suppliers to buy power from, 

commitment problem". 25 Note that suppliers who submit a bid with a high price may 

To maintain high reliability in the electric system, a utility end U P ™t selling any of their production. On the other hand, 
refines its load forecast every eight hours and solves, using selling power at a low price may not create enough revenue 
the updated forecast, a new unit commitment to adjust the f° r a generator. Given that no one knows in advance the 
schedule. The new schedule must respect all constraints amount of power that competitors may bid for, the electric- 
imposed by the previous one. For example, if a unit, that has 30 power market will become more uncertain and risky, 
a minimum up-time of eight hours, was switched on two The hope is that deregulation will result in cheaper prices 
hours before revising the schedule, the revision must make and improve the overall economy by encouraging invest- 
sure that this unit runs for at least another six hours even ments in electric utilities. The size of the electric industry is 
though this unit may not be needed. Clearly, such ineffi- expected to grow after deregulation as was the case with the 
ciencies may arise if the load forecast changes significantly 35 telecommunications industry. In the case of 
due to unexpected elements such as a sudden change in telecommunications, industry revenue shot up from $81 
weather patterns or generator failure. billion to $170 billion in ten years. 

For a more detailed discussion on the subject of power On the other end of the spectrum, due to deregulation, the 

generation and transmission, see A. J. Wood and B. E. 4Q load on a utility system is becoming increasingly unpredict- 

Wollenberg, Power Generation, Operation, and Control, able. The reason is trading transactions that change the load 

John Wiley, 1996. pattern significantly. Some utilities, for example, may sell 

The power industry is now going through deregulation. more than 30% of their power generation to other utilities on 

The current picture of a single utility controlling the market certain days. Demand and supply in the market are functions 

in a specific region will soon disappear. Instead, there will be 4S of volatile electricity prices which in turn depend on highly 

power producers who sell their production to a power pool; unpredictable elements such as weather conditions around 

and power suppliers who will buy power from the pool and the country and fuel prices which may vary within a wide 

sell it to their customers. Although the full picture of the range. 

power industry after deregulation is not yet known, it is clear TUUADV ™ TOC imraxmnM 

that utilities need to prepare themselves for an open market 50 SUMMARY OF THE INVENTION 

in which buying and selling power are to be considered It is therefore an object of the invention to provide a 

when the schedule of the generating units is created. computer risk-management system for scheduling the gen- 

The main reason behind deregulation is the high price of erating units of an electric utility, 

electric energy. The first step towards deregulation was It is another object of the invention to provide computer 

taken in 1978 with the passage of the Public Utilities 55 implemented process that manages generating units of an 

Regulatory Policy Act. This act encouraged non-utility gen- electric utility which handles multiple fuel, fuel constraints, 

eration and required utilities to buy power from independent varying fuel prices, power trading, and load uncertainty, 

generators. The Energy Policy Act of 1992 took deregulation The invention provides a computer implemented process 

a step further by mandating open access on the transmission for scheduling the generating units of a utility while taking 

system for wholesales. 60 into consideration power trading with other utilities and the 

Currently, electricity is sold as a service that is delivered stochastic load on the system. The system allows the user to 

at specified points. For example, each one of us expects to provide multiple load forecasts and to vary the fuel price 

receive electric power through a meter outside the house. We between the different scenarios and the different periods of 

also pay for this service regardless of its producer or which the planning horizon. The tool allows the user to model 

power lines it followed. That is, the electricity bill indicates 65 accurately the uncertain trading transactions and the chang- 

the total usage of electricity in KWH and the service price ing fuel prices. Given (1) a planning horizon, (2) a set of 

per KWH without incorporating any other details into the electric-load forecasts and fuel prices, (3) a full description 
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of the properties of the electric-power generators, (4) reserve taking power trading and fluctuation in fuel prices into 
requirements for the system, (5) an estimate of the price of consideration. This tool mixes traditional techniques used in 

electricity in the open market at each hour of the week, and scheduling generators of an electric utility with hedging 

(6) a set of possible trading transactions for the next two to strategies that are widely used in the finance industry. The 
seven days, the goal is to meet the electric demand of 5 result is a robust generating schedule, 

customers at a minimal cost while making the maximum FIG. 2 is a block diagram showing the relationship of the 

profit possible from power trading. To do so, a mathematical risk management system tool according to the invention and 

model of the problem is solved using appropriate optimiza- the generators of an electric utility. The tool 111 is part of a 

tion techniques. The solution provides the status of each control center 11 and receives as inputs fuel price and load 

generator at each time period of the planning horizon under 10 forecasts 112, projected trades 113, and generator data 114. 

each given scenario. By "status of a generator", what is The output of the tool 111 is the fuel consumption, fuel mix, 

meant is whether it is on or off. The solution also provides and generation requirement to each of a plurality of gener- 

the load on each generator during each period in which it is a ting plants 12j to 12„ which are connected to transmission 

operating, an optimal fuel mix for each generating unit, and ii nes 13 f or delivery of electrical power to customers 14 2 to 
the prices for purchasing and selling power in the periods of is I4 m and to other utilities 15. The inputs and output of the 

the planning horizon. The technique used to solve the model risk management system 111 are illustrated in FIG. 3 in the 

provides information regarding the sensitivity of the solution form of spreadsheets. 

to the input parameters and other valuable information to the Here there is introduced the notion of n0 n-anticipativity 

decision maker. because of its frequent use in the following sections. The 

BRIEF DESCRIPTION OF THE DRAWINGS 2 ° model lnat ^ used in tne following description is a multi- 
stage mathematical program with varying right-hand side. 

The foregoing and other objects, aspects and advantages By "multi-stage", what is meant is a model in which 

will be better understood from the following detailed decisions are made sequentially over a certain planning 

description of a preferred embodiment of the invention with horizon. The "varying right-hand side" reflects uncertainty 

reference to the drawings, in which: 25 m some of me svstem parameters. In the case of power 

FIGS. 1A and IB are graphs illustrating the relationship generation, the electric load on the system (which appears in 

between heat input and electric power output; the right-hand side of some of the constraints) is unknown 

FIG, 2 is a block diagram showing the relationship * n advance, 

between the risk management system according to the 30 0 Qe ca ° approximate the uncertainty in the right-hand 

invention and the generators of a utility; side using a set of scenarios each of which has a probability 

FIG. 3 is an illustration of input and output of the risk that reflects its importance or likelihood. The stochastic 

management system according to a preferred embodiment of P ro g ra m becomes that of minimizing the expected cost of 

the invention- operating the electric system under the given scenarios. A 

FIG. 4 is 'a flow diagram of a stochastic model for 35 J*"* 1 ? SOlu ;i°° for th * s '? chas,i * P ro f a ?? ^^Y^ 

determining generation levels under different conditions; Z Y ?K t° ^ N° D -anUcipa mty 

r-ii- e • a states that if two different scenarios, s and s f , are mdistin- 

FIG. 5 is a flow diagram of another stochastic model for guishable at time t on the basis of information available 

determining generation levels under different conditioas; about mem at ^ t> then the dec isions made for s and s' 

FIGS. 6A and 6B are graphs showing examples of a must be the same, 

heat-input versus power-output curve and the corresponding 40 To further clarify the concept of scenario bundles and 

generator cost function; their impact on decisions, an example from the finance 

FIG. 7 is a flow diagram of the overall operation of the industry is used. Assume that we want to invest $1,000 so 

risk management system according to the invention; that the return after three years is maximized. To simplify the 

FIG. 8 is a time line illustrating on-peak and off-peak time 45 mode l> we assume that there are two possible investments, 

periods; stocks and bonds, and that we can reallocate money from 

FIGS. 9A and 9B time lines showing the use of trade one investment to another by the end of each year. One way 

forecasts to create new scenarios* 10 model lms P rot) l em is bv using scenarios to represent the 

FIG. 10 is a flow diagram of the BUNDLE procedure ^lure retUmS °? ™* is, we study the market 

implemented in a preferred embodiment of the invention; 50 c ° me ^ r f™ rates for b ° th ^vestments 

. .. .. . (stocks and bonds) over the three years of the planning 

FIG. 11 is a time line showing a scenario tree with four horizon. The scenarios are described in Table 1. 
scenarios; 

FIG. 12 is a flow diagram of the SUCP procedure imple- TABLE 1 
mented by the preferred embodiment of the invention; 

FIGS. 13A and 13B are, respectively, a time line and a 
diagram showing the process of solving a dynamic program 
for a scenario tree with four scenarios and six time periods; 

FIG. 14 is a flow diagram of the DP procedure imple- 
mented by the preferred embodiment of the invention; and 60 

FIG. 15 is a flow diagram of the FUEL procedure imple- 
mented by the preferred embodiment of the invention. 
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6% 
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12% 
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DETAILED DESCRIPTION OF A PREFERRED ~ . . t . „ # 

CMnnnTKyrcKrr nn tue rxn/cMnnM 10 maximize the expected return, our investor must allocate 
EMBODIMENT OF THE INVENTION „ C1 , u , n r 4l jl,^, , 

65 the $1,000 between slocks and bonds. To do so, he or she 

The present invention is a tool that can be used to uses an (stochastic) optimization to develop a decision tree 

schedule the generating units of an electric utility while for the future. The decision tree indicates where to invest in 
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the beginning of the first year and how to reallocate the forecast. It is assumed that the load forecasts and their 
funds by the end of each year. The reallocation is a function weights are given in the form of a computer file, e.g., a 
of the return values observed during a year. For example, if spreadsheet, on the hard disk of a computer system, 
we observe a return of 10% on stocks and 5% on bonds in 2. The physical and operational properties (specifications) of 
the first year, then we know that Scenario 5 is out of the 5 the generating units: Each utility maintains accurate infor- 
picture and that the future decisions are based on scenario 1, mation regarding the specifications of its generators. The 
2, 3, and 4. At this point (end of year 1), the decisions made description of each generating unit can be summarized 
for scenarios 1, 2, 3, and 4 must be the same. The fact that using ten numbers: minimum operating capacity, maxi- 
we must adopt similar decisions for these four scenarios is mum operating capacity, minimum up-time, minimum 
called "non-anticipativity" and the corresponding decisions 10 down-time, start-up heat requirement, shut-down cost, 
are said to be "bundled". Let us take this example one step current status, and the heat-input versus power-output 
further and assume that we observe return rates of 8% and curve which requires three parameters. Again, it is 
7% in the second year. At the end of year 2, we know that assumed that this information is provided on the hard disk 
the only two possible future scenarios are Scenario 1 and in the form of a table or a spreadsheet. 
Scenario 2. The decision tree provides an investment deci- is 3. Fuel cost and its availability: To translate the curves of 
sion for these two scenarios that is independent of the FIGS. 1A and IB into cost, it is assumed the knowledge 
outcome of Year 3. Here again, we say that scenarios 1 and of the cost per BTU for each fuel used in generation (see 
2 are bundled at the end of Year 2. FIGS. lAand IB and FIGS. 6Aand 6B). This information 
The process described is implemented on a computer is needed to compute the start-up cost and the operating 
system for scheduling electric-power generators of a utility 20 cost of a generator. Our system permits the user to specify 
company. Given (1) a set of load forecasts on the electric different cost/BTU for a fuel for different hours over the 
system, (2) a prediction of trading transactions that may take course of the week and for different load scenarios. This 
place within the next 168 hours, (3) a description of the allows the user to change the prices of fuel to reflect his 
physical properties of the generating units and the expected or her view of the spot-market. The user can also provide 
spot market prices for electric power, and (4) fuel prices 25 upper and lower bounds (and any other linear constraints) 
associated with load forecasts, our process provides for each on the amounts of fuel that can be used at each hour or 
scenario (1) generation levels, (2) fuel usage, and (3) prices day. It is assumed that fuel information is stored on a 
at which power can be traded. This information is then used computer disk and is accessible by our system. The linear 
by decision makers and plant managers to run the electric constraints are also assumed to be stored on the disk in a 
system reliably. Our tool provides multiple decisions, e.g., 30 tabular form or in the standard Mathematical Program- 
generation, fuel consumption, and power prices, for each ming System (MPS) format described in IBM document 
hour of the planning horizon. Each of these decisions is number 

optimal under certain circumstances. The system adminis- 4. The projected power-trading transactions in the market: A 

trator observes the market, then looks up the appropriate view of the market and what type of trading transactions 

decision from a table provided by our tool. 35 may take place must be supplied by the user. For example, 

What distinguishes the system of this invention from the user may indicate that there might be an offer to "buy 

existing ones is the incorporation of fuel prices into the 200 MWH at S20/MWH tomorrow on-peak with a prob- 

model and allowing these prices to vary. The invention also ability of 10%" or a market to "sell 100 MWH at 

allows the user to impose constraints on fuel consumption. $25/MWH in two days off-peak with a probability of 

Incorporating fuel into the model requires iterating between 40 25%." It is assumed that this information will be supplied 

a stochastic unit-commitment solver and a stochastic linear- in a computer file that is stored on the hard disk, 

programming model. The stochastic linear program alio- Given access to the previous data, our tool produces a 

cates fuel optimally between the different generators of the table (spreadsheet) that indicates the load (MWH) on each 

system. The cost function of each unit is then constructed unit at each hour of the week. This table is passed on to the 

and a stochastic unit commitment is solved. The process is 45 administrators of each generator to be used as an operating 

repeated until a convergence is reached. schedule for the coming week. It also provides prices at 

The present invention provides a utility with a tool that which power is expected to be traded. As a result of the 

promotes a better understanding of the relationship between calculations, each plant manager is given precisely what 

the electric-power and fuel markets allowing a utility to type of fuel to use and how much to burn at each hour of the 

hedge against uncertainty in both markets. This is done by 50 week. The tool can be used to perform a "what-if * analysis 

modeling the customers' varying demand, the uncertain which helps the user in developing a better understanding of 

electric-power market, the generating units, the availability the relationship between the power market and his/her 

of different fuel, and the fluctuating fuel price in the open electric system; and in better quantifying the risk involved in 

commodity exchange. buying and selling power in the periods to come. 

Here is a more accurate description of the input of this 55 The existing tools for scheduling electric -power systems 

process. use a single load forecast when searching for optimal 

1. The forecasted electric load on the system of the utility: decisions. Furthermore, they do not consider power trading 

This is the load generated by customers who fall within in the model. In other words, existing tools assume that 

the utility's territory. Given load patterns observed in future loads are known in advance which results in conser- 

previous years and given load values for the past week, 60 vative schedules. What distinguishes our tool is that it allows 

utilities use the weather forecast (temperature, humidity, the user to incorporate risk, through predictions of the load 

sky conditions, and wind) to predict the load for the next and fuel prices, and uses these predictions to create optimal 

week or 168 hours. All utilities have advanced tools and schedules. Our tool uses hedging strategies to produce 

employ experienced planners to perform the forecasting robust schedules that minimize cost and manage risk effi- 

task. Due to uncertainty in load, the user is allowed to 65 ciently. 

provide a set of forecasts to better approximate the load The idea of using multiple scenarios to approximate 

distribution. The user also assigns a weight to each uncertainties in the load forecast was suggested in A. 
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Carpentier, G. Cohen, J. C Culioli, and A. Renaud, "Sto- 
chastic optimization of unit commitment: A new decompo- 
sition framework," IEEE Transactions on Power Systems, 
11(2): 1067-1073, May 1996, and in Takriti, Birge, and 
Long, "A stochastic model for the unit commitment 
problem," IEEE Transactions on Power Systems, 11(3): 
1497-1508, August 1996. However, both models assume 
that fuel prices are known in advance and are fixed over the 
planning horizon. 

FIG. 4 is a flow diagram illustrating the stochastic model 
of Carpentier et al. Load forecasts and generator properties 
are input in function block 41. Then, using Lagrangian 
relaxation, the problem is decomposed into single-generator 
subproblems in function block 42. Lagrange multipliers are 
initialized in function block 43 prior to entering a processing 
loop. In the processing loop, the stochastic dynamic program 
of each generator is solved in function block 44 and, after 
each solution, a test is made in decision block 45 to 
determine if the current solution is optimal. If not, the 
Lagrange multipliers are updated in function block 46, and 
the process loops back to function block 44. When an 
optimal solution is obtained, the generation levels under 
different scenarios are output at function block 47. Carpen- 
tier et al. do not discuss issues related to the practicality of 
their model and do not report computing requirements and 
execution-time for their technique. 

The flow diagram for the stochastic model of Takriti et al. 
is shown in FIG. 5. The load forecasts and generator 
properties are input in function block 50. Lagrangian relax- 
ation is used in function block 51 to relax the bundle of 
constraints. In this case, each scenario has its own problem. 
Progressive hedging penalties are initialized in function 
block 52 before entering the processing loop. In the pro- 
cessing loop, Lagrangian relaxation is used in function block 
53 to decompose the problem into single-generator subprob- 
lems. Then, a subprocessing loop is entered beginning with 
initialization of the Lagrange multipliers in function block 
54. The deterministic dynamic program from each generator 
is solved in function block and, after each solution, a test is 
made in decision block 56 to determine if the current 
solution is optimal. If not, the Lagrange multipliers are 
updated in function block 57 and the process loops back to 
function block 55. When the current solution is optimal as 
determined in function block 56, a further test is made in 
decision block 58 to determine if bundle constraints are 
satisfied. If not, the progressive hedging penalties are 
updated in function block 59, and the process loops back to 
function block 53. When the bundle constraints are satisfied, 
the generation levels under different scenarios are output at 
function block 60. Takriti et al. use a parallel machine with 
thirty-two nodes to handle the large size of their model. 
Given that a comprehensive model has a significantly larger 
size than the examples of Takriti et al. and that a parallel 
platform is very costly, this technique may be hard to apply. 

The system according to the present invention builds on 
Carpentier et al. and Takriti et al. but uses different math- 
ematical techniques to accelerate the process of solving the 
model and to make it more general. To be more specific, the 
present invention implements: 

1. A storage scheme for the scenario tree and bundle 
constraints: The result is significant reduction in calcula- 
tions and storage requirements. 

2. A new updating mechanism for the Lagrange multipliers 
of the stochastic unit commitment problem: the bundle 
trust-region method is used. 

3. A clear and concise treatment of power trading. 

4. A flexible operating range: the operating range is allowed 
to be a function of time which allows incorporating 
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maintenance into the unit commitment. Generators are 
permitted to have a varying capacity when they are turned 
on (ramping). 

The resulting system is a practical tool that runs efficiently 
on any computer platform and under any operating system, 
including parallel computers. Parallelism reduces the sys- 
tem's execution time (approximately) proportionally to the 
number of available processors. Given that the market 
profile may change quickly, solving the model in parallel 
permits the system administrator to respond efficiently to the 
dynamic market. 

System's Input 

First, there is discussed the input parameters required for 
the model of the present invention. The duration of the 
planning horizon is how far ahead we want to look so that 
our present decisions are close to optimal. A planning 
horizon varies from three to seven days depending on the 
utility. The planning horizon is divided into periods. Usually, 
each period is one hour. Our model assumes that the number 
of periods, T, is supplied. 

Second, planners must provide load forecasts (MWH) for 
each period of the planning horizon for the region covered 
by the Utility. Presently, all utilities employ trained indi- 
viduals to produce reasonable forecasts. There are also 
special forecasting packages that are widely used by utili- 
ties. The forecast is a function of the weather in the region 
of interest, time of day, day of year, and many other factors. 
Methods for developing load forecasts are not discussed 
here. It is assumed that a set of S load forecasts is provided. 

Each forecast (scenario) provides the total load, d/, on the 
system at each time period, t=l, . . . , T, of the planning 
horizon. The superscript s indicates that a variable or a 
parameter is associated with forecast s. For example, d 2 4 is 
the total load on the system at time period 2 under scenario 
4. 

A weight, p f , is associated with each scenario, s, to 
approximate the likelihood of its occurrence in the future. If 
p^p^., then we are more likely to observe scenario s than to 
observe scenario s\ If p^O, then scenario s cannot occur, 
and we can drop it out of our scenario set. Note that 
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represents the probability that scenario s occurs. Table 2 
gives an example of load forecasts. 

TABLE 2 

Load on our Utilily, d/ 





P. 


t - 1 


2 


3 


4 


5 


s - 1 


0.35 


10495 


11135 


11653 


11829 


11995 


s - 2 


0.56 


10846 


10982 


11563 


11674 


12039 


s - 3 


0.09 


11138 


11956 


12310 


12388 


12583 



To reflect uncertainty in fuel prices, the user is allowed to 
provide different prices for a fuel, j, over the different points 
in time, t, and under the different load scenarios, s. The price 
of fuel j is denoted by 
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quadratic cost functions of the individual fuels. Given that 

/ $ \ we switch from a cheap fuel to a more expensive one as we 

'MbtuJ' increase the load, the function i t ' is convex which makes 

the task of minimizing it a simple one. From here on, it is 
5 assumed that the cost function, f ( / (x) is convex and 

Given that the goal is to minimize the expected cost of increasing and that f ( . / (g i>t )>0. It is our experience that cost 

running the electric system, the expected value of fuel cost functions faced in power generation fall within this class of 

over scenarios that belong to the same bundle can be used. functions. Hence, our assumption does not pose any limi- 

That is, if Q is the set of scenarios that satisfy the condition tations on the applicability of this technique. It is also 

d tt -d„. for any s and s' in Q and for t-1, . . . , t, then the assumed that we have a tool for finding the minimum of 

expected value E jen (h,/)can be used as the cost of fuel j at 10 f. /( X( . f)'X t %/ over the range [g, jf G,-,] where X/ is a 

time t»l, . . . , x under the scenarios of Q. non-negative real number. Since f lV J 'is convex, the previous 

Third, the model requires a complete description of the minimization is trivial, 

generating units of the utility. By a "complete description", (c) Minimum Up-time and Minimum Down- time: When a 

what is meant is providing the total number, N, of generators generator is switched on, it has to remain on for at least 

in our system and the operational properties of each gen- 15 L _ p^ods and when it is switched off, it has to remain off 

erator. That is, for each generator, i e {1, . . . ,N}, the user for at j east 1( . periods. Note that if either of these two 

must provide: constraints does not exist, one can set the corresponding 

(a) Operational Range: This is represented by the minimum parameter to zero. For example, if there is not a minimum 
power, fcfff and the maximum power, G S/t that can be up _ time requirement on a generator, i, then L, can be set 
produced by this generator during penod t. Clearly, 20 t0 0 and the C0DStraint ^ dropped. The values of L, and L 
and must satisfy the constramt O^g, S G „. When it is must be non . negative intcger for & { e{1> . . . >N} . 

on (during period t), a generator, i, can produce any (d) Start and shut-down Costs: A generator must be 

amount of power in the range [g,„ G f J. Note that the warmed fof a number of periods before it starts 

operational range depends on the time t. Although, the generating electricity. The cost of fuel (usually, natural 

minimum and maximum operating capacity for a genera- 25 gas) ^ at this stage ^ called a start . up cost It ^ a 

tor are not affected by time, t is incorporated to allow Qf the faeatj gu * BTU( r ifed tQ wafm ^ 

considering maintenance in the model. For example, if a UQ - t and tbe - ce 

unit i is scheduled for maintenance over a certain time 

range, then its capacity over this range can be reduced or $ 

made zero (it is not always necessary to shut down a unit 30 ^(btu)' 

during maintenance). 

In some cases, a unit cannot operate at its maximum 
capacity immediately after it is started. For example, when of the heating fuel during the warm-up period. The start-up 
the load reaches g i t (or some other level), it must be kept at cost is denoted by SU,/. A shut-down cost, SD,/, is 
that level for a few periods (usually two hours). Such limits 35 incorporated in the model. Most utilities set SD ( -/ to zero, 
are applied directly to the states of the dynamic program of (e) Status of Each Generator Prior to the Beginning of the 
each generating unit. Planning Horizon: Given the minimum up-time and mini- 
Note that we can allow g, , and G i t to vary under different mum down-time constraints, we need to know for how 
scenarios in the model. However, there is no need for such long a generator has been up or down prior to the first 
flexibility in practice. 40 period. The reason for that is to make sure that we satisfy 

(b) Heat-input versus Power-output Curve: This is the curve minimum up-time and down-time constraints at the begin- 
that represents the relationship between the heat in the ning of the planning horizon. 

burner and electric-power output (see FIGS. 1A and IB). (f) Classification of Generators: Generating units are usually 
This curve is denoted by f/ 1 where i e {1, . . . , N}. In other divided into different classes. Some units are classified as 
words, to generate x units of power (MWH), the burning 45 "must -run units". As the name indicates, these units are 
fuel must provide f *(x) thermal units (BTU). The func- functional all the time and we do not need to worry about 
tion f* is quadratic, convex, and increasing. It can be fully scheduling them. An example of such units are nuclear 
described by providing three parameters: a f , b ( , and c t -. generating units. Another class is the "cycling units". 
That is, the heat required to produce x units of power for These units are switched on and off during the planning 
one period using generator i is a i -x 2 +b,-x+c t . To guarantee 50 horizon depending on the load on the system. A third class 
that f * is convex and increasing, the parameters a f , b £ and is "speakers" that are used during emergencies and high- 
er must satisfy a^O, -0.5 b/a^g^ if a t >0, b f -^0 if a/), and demand periods. Finally, we have the "unavailable" units. 
a *g *> +D t&> +c i=°* These are units scheduled for maintenance or unavailable 
The cost function of generating electricity on generator i for other reasons. Knowing the types of units helps 
is denoted by f ( -/. Although t/ 1 is a quadratic function, the 55 accelerating the solution procedure, 
function t ir s may not be quadratic. The reason is that when (g) Loss Factors: To consider losses in transmission lines, 
the upper capacity, G irL , is approached, a generator usually one can use a loss factor for each unit to adjust its 
requires a more combustible fuel. Some generators, for production. For example, if the average transmission 
example, require the use of natural gas in addition to coal losses for a unit is 9%, then one can apply a reduction 
when the load exceeds 80% of G ir As a result, the geoer- 60 factor of 91% to the production of this unit. Note that 
ating cost curve must combine the cost curves of both fuels. these factors can be incorporated simply by multiplying 
FIGS. 6 A and 6B show an example of a cost function f,-/. the heat curve of a unit by a constant that depends on its 

The generation levels at which we switch fuel is denoted penalty, 

by G tj , , k=l, . . . , K, where K+l is the total number of Fourth, to avoid any blackouts, utilities make sure that, at 

different fuels (excluding the start-up fuel) used on i. The 65 each period, the maximum operating capacity of their sys- 

values of JS t f are assumed to be known for each generator. tern exceeds the demand of this period by a certain amount. 

The cost function f ( / can then be constructed from the This excess capacity is called "spinning reserves". To clarify 
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the concept of spinning reserves, assume that our system has generator may use more than one fuel at any time period), 

ten generators of which seven are on line at the current time An easy way is to use the fuel with the least cost as an 

period. Let us also say that the forecasted demand at this estimation of E,- then update this value (as in function 

period is 12,000 MWH. In an optimal solution, the total block 78). In this way, the cost function, f,/, of generating 

generation would be close to 12,000 MWH. However, if the 5 a MWH of power on generator i is obtained by multiplying 

maximum capacity of our operating units; i.e., the total f* by E,/. 

capacity when each unit is operating at G i f , is close to A processing loop is now entered which begins by calling 
12,000 MWH, our system cannot take any unexpected the procedure SUCP in function block 75 to obtain the load 
increase in the demand. In other words, the reliability of our on each generator under any scenario. Given f lV % a stochas- 
system is low. The reliability can be improved by forcing the 10 tic unit-commitment model is solved. As a result, the gen- 
total maximum capacity of the operating units to exceed the eration of each unit at each time period (under the different 
expected demand by a certain amount of power. This excess scenarios) is obtained, and this can be used to compute the 
capacity is the spinning reserves and is indicated by r„. heat-input needed for each unit throughout the planning 
There are other reserve constraints that can be enforced to horizon. Then given the load on each generator, the proce- 
improve reliability. The treatment of such constraints is very is dure FUEL is called in function block 76 to determine an 
similar to our treatment of spinning reserves. optimal fuel allocation. 

Fifth, to incorporate trading transactions, the user of our To satisfy the required neat mput for each generator, a 

system must provide his or her view of the future power pr0 g ram i s solved in which the decision variables are 

market. Using the weather forecast and information that may the amount of fuel burnt on each unit at eacn time period of 

influence power usage in other regions of the country, power 20 course> our decision variable has a component for each 

traders estimate the price of buying or selling power in the possible fuel in the model. The linear program is solved 

market and the quantities that may be offered for trading. su5 j ect to (he constraint that fuel fed to a generator must be 

The projection of trading transactions is performed for the sufficient to generate the heat needed. It must also consider 

next two to seven days. The tool according to the invention non-anticipativity and other fuel constraints imposed by the 

uses four parameters to describe a transaction. These param- 25 user solution of the linear program provides the fuel 

eters are: amount of power, v(j) MWH; price per unit, % quantities that need to be fed to each generator at each hour 

S/MWH; the probability, that such a transaction may take of the pi ann ing horizon so that the (expected value of) cost 

place; and the delivery time of the transaction. The delivery ^ minimized. 

time of a transaction indicates the date and time at which * * * ■ * j • j * • ui i nn* j ♦ • r*t. 

•ii . j i- i ™ , i. . , , A test is next made m decision block 77 to determine if the 

power will be delivered. The delivery time can be on-peak 30 . A , • f . „ A _ . ^ _ Tf . n M ,„ 

K % „ ™ . „ ^ / «• i /. current fuel allocation is the same as before. If not, a new 

(between 7:00 a.m. and 10:59 p.m.) or off-peak (between * a *• r u * . j • c 

i„ ^ * ™* i * -j 1 «• cost function for each generator is computed in function 

11:00 p.m. and 6:59 a.m.). Weekends are considered off- , t , - 0 . r ° , . , 4 £ , , , 

, y ' block 78 before the process loops back to function block 75. 

P e t/ i -1- . • * i- . . Using the linear program solution in function block 78, a 

Currently, a utility uses the information outlined above as r * r u . u * j u 

.-j j * i • jj- -ii j r new h it for each generator can now be computed by 

the first, third and fourth in addition to a single load-forecast 35 avera £ me values of 
from the second to generate a reasonable schedule for its 
generating units. In many cases, the effect of fuel on a 

generator's cost function is approximated and fuel con- h*{ ] 

straints are ignored. The multiple forecasts of the second and IBTUJ 

the uncertain nature of trading transactions and fuel prices 40 

are yet not incorporated into generation in the power indus- over all fuels fed to a given generator. The new value of E f / 

try- yields a new cost function, f,- /, that is used in resolving the 

Solution Technique unit commitment problem ' 

As the process is repeated, the fuel mixtures and the load 
The overall structure of the risk management system 45 0 n the generators approach an optimal solution for the 
according to the present invention is shown in the flow problem. If the test in decision block 77 is yes; i.e., the fuel 
diagram of FIG. 7. Load forecast, trading forecast and allocation has not changed as a result of the new iteration, 
generator and fuel properties are input in function block 71. then the process moves to function block 79. In this case, the 
The first step in processing the input data is to generate load process has converged to an optimal solution. Function 
scenarios. These scenarios are a result of applying power- 50 block 79 outputs the results of the final iteration. This 
trading contracts to the original load forecasts. The resulting includes the status of each generating unit at each time 
scenario tree approximates future uncertainties in load and p er j od and under each scenario, the amount of power 
trades. The probability of each scenario is computed by generated on each generator at each time period and under 
multiplying the probabilities of the trades with those of the eacn scen ario, the marginal prices or the Lagrange 
load forecasts. Scenarios with relatively small probabilities 55 multipliers, and the fuel usage provided by the linear pro- 
are eliminated and their probabilities are distributed among gram solution. The process is then terminated, 
other scenarios. The scenarios are generated by calling the 

procedure TREE in function block 72. Procedure TREE is Procedure TREE 

described in more detail below. In some cases, the user may The procedure TREE is called in function block 72 in the 

provide a scenario tree which must be processed to define 60 process shown in FIG. 7. Given a forecast for power-trading 

the bundle constraints. Recall that two scenarios are mem- transactions in the next two to seven days, TREE generates 

bers of the same bundle at time t if their loads are the same a scenario tree that represents trade opportunities. Note that 

up until t. Thus, if needed, the procedure BUNDLE is called we do not concern ourselves with trading transactions on the 

in function block 73. following day because utilities are required to close all deals 

Using the fuel available for consumption by each 65 twenty-four hours in advance, 

generator, i, we compute, in function block 74, an approxi- To create a scenario tree, the user supplies his or her 

mate value for the cost per BTU or H,. / (recall that a vision of the power market in the next two to seven days. A 
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"vision" is represented by the size of a contract in MWH, the 
cost of contractual power per MWH, the type of contract; 
i.e., whether it is buying or selling, the day on which the 
delivery will take place, and the time of delivery. The time 
of delivery can be cither on-peak or off-peak. The on-peak 
period is between 7:00 a.m. and 10:59 p.m. while the 
off-peak period covers the rest of the day. Weekends are 
considered off-peak. The user must provide a probability for 
each contract. FIG. 8 is a time line illustrating on-peak and 
off-peak time periods. 

Given that we may initiate a trade at two points in time for 
a given day, the scenario tree branches at most at two points 
on all days except for the weekend. On Saturdays and 
Sundays we have one branching point per day. To create the 
scenario tree, it is assumed that the contractual load for the 
first day is known (because the contracts are closed for this 
day). For day two, a number of branches are created equal 
to the number of off-peak contracts. The branching point is 
the beginning of day two. The probability of each branch is 
the same as the probability of the corresponding trade. At 
7:00 a.m. of day two, for each branch of the tree a set of 
branches are created to reflect the on-peak potential con- 
tracts. The probability of each new branch is equal to the 
probability of its contract times the probability of its pre- 
decessor. If the probability of a branch is below a certain 
threshold, e, the branch is eliminated. The probability of this 
branch is distributed among other branches that share the 
same predecessor in proportion to their probabilities. The 
process is repeated until the last branching point in the 
system is reached. 

The previous process may yield a large number of sce- 
narios. In practice however, due to the difficulty in predict- 
ing the market on a horizon longer than three days, power- 
trading forecasts are given for day two and day three only. 
Another factor that reduces the number of scenarios is the 
use of e. A scenario with a probability of e or less is assumed 
to be insignificant. Our experience indicates that in a rea- 
sonable tree, e can be set to 2% without affecting the hedging 
results. In other words, it is sufficient to consider fifty or less 
scenarios. 

The user may prefer to provide a set of load forecasts (and 
a set of corresponding fuel prices h /J5 ) instead of a single 
one. In this case, the process of branching, using trade 
forecasts, needs to be performed on each load scenario. The 
probabilities need also to be multiplied by the probability of 
each load forecast. To give an example, consider the sce- 
nario tree of FIG. 9A. The user provides two possible 
forecasts. The first has a probability of x h and the second has 
a probability of Yz. Let us say that we expect to sell a certain 
amount of power with a probability l A on the beginning of 
period four. Then, each load forecast branches at the begin- 
ning of period four creating two branches. The final prob- 
abilities are shown in FIG. 9B. 

A depth-first search is used to generate the scenario tree. 
The steps required when one load forecast is supplied are 
described. Hie method can be generalized easily to the case 
in which multiple load forecasts are provided. To fix 
notations, assume that the load forecast is H(t), t-1, . . . , T. 
Denote the number of branching points by M and the time 
periods at which a branching takes place by x(i), i*»l, . . . , 
M. Q(i) is defined to be the set of all trades that can start at 
time T(i). The user of our system supplies for each trade, 
j=l, . . . , J, its end point, xfy); its volume, v(j) MWH; cost 
per unit, £(j) $/MWH; and its probability, Jt(j). Note that v(j) 
is negative if j represents a purchase offer and positive if j 
is a selling offer. 
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The approach is clarified through the example of Table 3. 



TABLE 3 
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71 
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0.50 
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20 


25 


0.25 


3 
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10 


15 


0.25 
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0.67 
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4 


-10 


10 


0.33 



Assume that we have four periods. We can sell twenty MWH 
with probability Ya or sell ten MWH with probability Va. 
Each of these two selling contracts starts at the beginning of 
period two and has a duration of two periods. At the 
beginning of period four, we may buy ten MWH with a 
probability of VS. Note that Q(1)={1,2}, Q(2>{3}, t(1>2, 

and x(2)~4. It is assumed that 3(t)=0, t-1 4. Depth-first 

produces the scenarios as shown in Table 4. 



TABLE 4 



s/t 
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Az 
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0.33 


1 
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4 
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-10 


0.17 


4 
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100 
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20 


20 
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0.17 
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-500 
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20 


20 


-10 


0.08 
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-400 
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10 


10 
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0.17 
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12 


-150 
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10 


10 


-10 


0.08 
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13 


-50 



To manage the scenarios efficiently, for each scenario, s, the 
first point in time at which it branches off the tree, x(s), its 
predecessor, o(s), and the change in the objective function 
cost due to buying/selling power, Az(s), are stored. The 
scenarios are stored in a string of real numbers, d. For 
example, the scenarios of Table 4 are stored in the form 
(0,0,0,0,-10,20,20,0,-10,10,10,0,-10). The variable k(s) 
indicates the location in which the last element of scenario 
s is stored. For example, Scenario 1 is stored in d(l) . . . 
d(ic(l)). For Scenario 2, only the fourth period needs to be 
changed. Hence, its elements are stored in d(ic(l)+l . . . k(2)) 
or 6(5). This scheme reduces the memory needed to store the 
scenario tree. If we guarantee that there is no overlap 
between the trades of two different branching points, it is 
possible to further refine the previous approach. 

The following is a pseudo-code of the depth-first search. 
Input, integer T, integer M, integer J, float 3(1 . . . T), integer 

x(l. . . M), integer Q(l . . . M,l . . . |Q(i)|), integer x^l . . . 

J), float v(l . . . J), float £(1 . . . J), float rc(l . . . J) 
Output, integer S, integer x(l . . . S), integer a(l . . . S), 

integer k(1 . . . S), float Az(l . . . S), float d(l . . . TxS) 
Local Variables, integer i, integer k, integer 1(1 .. . M) 
Initialization. 

1. Set k«-M, S*-l, and 1(0—1, ■ ■ . M- 

2. Set x(S)«-l, o(S)«-l, and k(S)«-T. 

3. Add v(Q(iJ(i))) to H between x(i) and X/(Q(i,l(i))) for all 
i-1, . . . ,M. Store the resulting vector in d(l . . . k(S)). 



60 



4. Set £lz(S) «- £ v(£l(l, KODflW /(/))) 



General Step. 

1. Increment l(k) by 1. 
65 2. If l(k)=|Q(k)|+l, then 

set l(k)«-l and k+-k-l. 
If k-0, then terminate. 



06/22/2004, EAST Version: 1.4.1 



6,021,402 



17 



18 



Go to 1. 

. For each i=k, . . . , M, add v(Q(l(i))) to 3 between x(i) 
and x,((fl(l(i))). Store the resulting vector in 
d(K(S)+l . . . K(S)+l+T-x(k)). 

. Se^a(S+l)«-S, x(S+l)«-x(k), and k(S+1)«-ic(S)+1+ 

T-x(k). 



5. set + 1) «- £ vwamwm- 



6. Set S*-S+l and k«-M. 

7. Go to 1. 

At termination, one can eliminate all scenarios that have 
probabilities less than e and distribute these probabilities 
among other scenarios. 

Procedure BUNDLE 

The procedure BUNDLE is called in function block 72, if 
needed, in the process shown in FIG. 7. The flow diagram of 
BUNDLE is shown in FIG. 10. Procedure BUNDLE begins 
in function block 101 in which the predecessor of Scenario 
1 is set to 1, its split time is set to 1, and the scenario index, 
s, is set to 2. In function block 102, the load of scenario s is 
compared to the load of scenario 1 to determine the split 
time of s. That is, we find the first point in time at which the 
load of Scenario s is different from that of Scenario 1. At this 
point, the predecessor of s is set to 1. A new variable, s', is 
introduced and initialized to 2. The idea is to compare the 
load of s to all other scenarios s'-l, . . . , s-1, in order to 
determine the split time and predecessor of s. The compari- 
son between Scenario s and Scenario 1 is performed in 
function block 102 while the comparison between s and each 
s' in {2, . . . s-1} is performed in function block 105. 

In decision block 103, a determination is made as to 
whether s' is less than s. If s'.iLs, i.e., s'=s, then the execution 
is transferred to function block 107 in which s is incre- 
mented by 1; otherwise, the execution moves to decision 
block 104. 

In decision block 104, the split times of s and s* are 
compared. If they are not equal, then the execution is 
transferred to function block 106. The fact that the split 
times are not equal indicates that s' should not be used as a 
predecessor to s. If the split times are equal, then s* might be 
used as a predecessor to s. In this case, the execution is 
transferred to function block 105. 

In function block 105, the loads of s and s' are compared 
to obtain a new split time for s. As one can see, the split time 
of Scenario s is getting higher as we go through more 
scenarios s' which results in lesser storage for the scenario 
tree. The predecessor of Scenario s is set to s\ Then, in 
function block 106, the value of s' is incremented by 1. This 
is how all scenarios from 2 to s-1 are scanned to determine 
the predecessor and split time of s. In function block 107, the 
index s is incremented by 1, as already mentioned. This is 
how all the scenarios in the scenario tree are scanned. 
Finally, in decision block 108, s is checked against the total 
number of scenarios, S, in our model. If s=S, the execution 
is terminated; otherwise, the execution is transferred to 
decision block 102 to continue processing. 

Thus, given a set of load forecasts, d/, t«l, . . . , T, 
s«*l, . . . , S, the procedure BUNDLE organizes these loads 
in the form of a tree. For example, if S«2 and T«3 and if 
d*-(l, 2, 2) and d 2 «(l, 2, 3), then d 1 and d 2 form a tree that 
has one branch at t-1 and t-2 and two branches at t-3. The 
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role of BUNDLE is to study the forecasts, d\ . . . , d*, and 
to define the different branches of the tree. Note that if the 
scenarios are generated by TREE, then there is no need for 
using BUNDLE. The procedure BUNDLE is provided to 
handle the case in which scenarios are provided by the user. 

To clarify the previous points, the scenario tree of FIG. 11 
is used. The process is initialized by setting x(l)«l and 
a(l)»l. Scenario 2 has the same demand as Scenario 1 for 
the first two periods. Then, x(2)=3 and a(2)=l. To place 
Scenario 3 in the tree, Scenario 3 and Scenario 1 are first 
compared which results in x(3)=»3 and o(3)=l. Now, given 
that x(3)«x(2)=3, Scenario 3 and Scenario 2 are compared 
starting at period 3. Scenario 3 branches off Scenario 2 in the 
fourth period. That is, we set x(3)=4 and o(3)=2. Finally, 
comparing Scenario 4 with Scenario 1 yields x(4)=2 and 
a(4)»l. Given that x(4)*x(2), we do not compare scenarios 
2 and 4. In the case of scenarios 3 and 4, we have x(3)*(4) 
and the process terminates. 

The following is a pseudo code of the procedure 
BUNDLE. 

Parameters, float e-10.0 

Input, integer T, integer S, float 3(1 ... S, 1 ... 1} 
Output, integer x(l . . . S), integer o(l . . . S), integer k(1 . . . 

S), float d(l . . . TxS) 
Local Variables, integer t, integer s, integer s\ 
Initialization. Set x(l)*-l, o(l)«-l, and k(1)«-T. 
General Step. For each scenario, s=2, . . . , S, 

Find the first time-period, t, at which |3(s, t)-d(l, t)|>e. Set 

t(s)*-1 and o(s)*-l. 
For each scenario s'=*2, . . . , s-1, if x(s)^r(s'), then find the 
first time-period, t>x(s), at which |d(s, t)-d(s', t)|>e and 
set x(s)*-t and a(s)«-s\ 
Set k(s)«-k(s-1)+T-x(s)+1. 

Procedure SUCP 

The procedure SUCP is called in function block 75 of the 
flow diagram shown in FIG. 7. The flow diagram for the 
procedure SUCP is shown in FIG. 12. The execution starts 
at function block 121 in which the Lagrange multipliers are 
initialized. There are different methods for this initialization 
in the literature. The following simple procedure is sug- 
gested: 

1. Using an average fuel cost, 



construct cost functions f, /. 

2. For each scenario and period in the planning horizon, 
Rank the generators in the ascending order of average 
marginal cost (S/MWH) at the maximum production 
level, G ( . f . 

Select enough generators so that the load and spinning- 
reserve constraints are met. 

Solve an economic power dispatch. The dual variables, 
X, and //, can be used to start the bundle trust-region 
method. 

Refer to S. Takriti, J. R. Birge and E. Long, ibid., for more 
details. 

In function block 122, the bundle-trust method for solving 
the dual problem is used. This method is described in H. 
Schramm and J. Zbwe, "A Version of the Bundle Idea for 
Minimizing a Nonsmooth Function: Conceptual Idea, Con- 
vergence Analysis, Numerical Results", SIAM Journal on 
Optimization, 2(1), pp. 121-152, February 1992. A brief 
description of this algorithm is given below: 
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Let n^)^^), where X 3 and fa are the dual multipliers 1. Solve a minimization problem for each period of the 

obtained from the Initialization step. Denote the Lagragian planning horizon. Each problem has two constraints: 

by £(jr) and a subgradient for it at k x by <x> v Set k«-l. demand and spinning reserves. These two constraints 

1. Solve the program m&xz{min 1SiSk <n i 7 ii+ti)?'(]z r K i d are relaxed using X and pi and the Lagrangian is 

5 optimized over the resulting period. 

t r, , r, v, 1 =r= 2. If the total generation cost is equal (close enough) to 

max- {mill! siak oj/5f + a>/ + £(*;)} 4- — W H. f » * • * • t_ • 

* 2t k that of the previous iteration, then terminate. 

3. Using the load on each generator, call the procedure 

2. If £(n k +7z) is sufficiently larger than £(x^, then either 10 FUEL. 

(a) enlarge i k and go back to 1; or 4. Compute a new cost function. 

(b) set re^i-jvwt, 5. Go to 1. 

(c) compute a new subgradient (o^ by calling DP for all p or a more detailed description of the economic power 
generators and computing the slack of demand and dispatch problem, the reader is referred to the book of A. J. 
reserve constraints, 15 Wood and B p Wollenberg, ibid. 

(d) increase k by 1 and go to 1. Summarizing, given a set of generators, {1, . . . , N}, and 

3. If £K+Jt) is not sufficiently larger than £(%), then a u re i a ted data, such as f,/, g,- „ G lV , and G i n and given a set 
either of i oa d forecasts on the electric system at each period in the 

(a) reduce t* and go back to 1; or planning horizon, SUCP solves a stochastic unit commit- 

(b) set Ttjt+^Jtjt, _ 2Q men t problem. That is, it produces, for each scenario, a 

(c) compute a subgradient, at %mi, schedule indicating the time periods in which each gener- 

(d) go to 1. The algorithm terminates when the norm of ami g unit is running and its load during these periods. The 
the subgradient is small enough. At this point, an resulting schedule guarantees that the total generation cost is 
optimal solution for the dual problem has been close to minimal and that the system provides enough supply 
found. Note that at termination of function block 25 t o meet demand. An important output of SUCP is the 
122, there may not be a feasible solution for the marginal cost, Xf ($/MWH), of power at each period. Under 
stochastic unit commitment problem. In function a given scenario, the marginal cost at time t is the cost, over 
block 123, the solution of function block 122 is used a ll operating generators, of producing one unit of power 
to search for a feasible solution for the primal during this period. Marginal costs are of great value to the 
problem; i.e., the stochastic unit commitment model. 3Q decision maker because they provide pricing information 
The preferred approach is that of F. Zhuang and F. D. that can be used in decisions related to buying and selling 
Galiana, "Towards a More Rigorous and Practical power. 

Unit Commitment by Lagrangian Relaxation", IEEE SUCP assumes that the cost curves, f,/, are given and 

Transactions on Power Systems, 3(2), pp. 763-773, independent of fuel prices (which are dealt with in FUEL). 

1998. The following is a brief description of their 35 Mathematically, SUCP minimizes the cost function 
approach: 

1. If all reserve constraints are satisfied, then proceed to s r (l) 

the economic power-dispatch phase. ^ + ^ ~ ^ l)+SUf ' ~ (w * " "W 51 * 



Weight the violations of reserve constraints by their 
probabilities. Find scenario s and period t for which the 40 
weighted reserve constraints are most severely vio- 
lated. wnere (u l -/-u (> _ 1 T (u t /-u tV _/)- are the positive and 
3. Increase by and resolve DP. negative parts of ^/-u^./. The minimization of (1) is 
4 Go to 1 subject to the demand constraints 

After obtaining a feasible solution, the execution moves to 45 
function block 124 in which an economic power dispatch is 
solved. This is the problem of allocating the electric load 
optimally between the different generators at each time 
period and under each scenario. Note that the economic 

dispatch problem assumes that units' status are given (from 50 and spinning-reserve constraints 
solving a unit commitment). Economic dispatch is a simple 
mathematical program in which each time-period is treated 
independently of others (unless we want to consider ramping 
rates). To solve such a problem efficiently, one can relax the 
demand and spinning-reserve constraints for each period and 55 

solve the dual. A good starting point is the Lagrange in addition to minimum up-time, minimum down-time, and 
multipliers obtained from the unit-commitment solver. operating capacity of each generator. Note that it is desirable 
Given that the generators' cost functions depend on fuel in some cases to impose additional constraints to improve 
prices, we revise the fuel mix using the economic power- the reliability of the solution. The treatment of any set of 
dispatch solution. If the changes in fuel mix are significant, 60 linear constraints that links the generators is similar to that 
we resolve the economic power dispatch model. Note that suggested for the above constraints, 
one can avoid iterating between economic power dispatch The procedure described is based on Lagrangian relax- 
and fuel assignment by merging- these two problems into a ation. To solve (1), the demand and spinning-reserve con- 
large-scale mathematical program. The standard economic straints are relaxed by associating non-negative multipliers, 
power dispatch method is modified to fit the model imple- 65 \ and //, with them. The value of X/ represents the marginal 
mented in the present invention. Here is the resulting pro- cost of producing power during t. Similarly, ft* represents 
cedure: the cost of having spinning reserves on-line at time period t. 



jTu?,G lV + = 1 T,s=l s, 
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The problem becomes that of maximizing the resulting must provide enough fuel to generate the required electric 

concave dual. This is the main attraction behind this power. That is, the amount of fuel used must satisfy the 

approach: substituting a difficult primal by a convex pro- constraint 

gram. On the other hand, the resulting Lagrangian is not 

smooth. The bundle trust-region method of Schramm and 5 ^ 

Zowe, ibid., is used to solve the dual. >?u = ' = 1 Nt 1 ~ 1 T> s = 1 s * 

Relaxing the demand and spinning-reserve constraints J " 1 
decouples the generating units. Each generator's problem is 
solved using dynamic programming. The solutions to the . 

dynamic programs provide the load on each generator. 10 where v ^.=° 1S the heat generated (BTU) by fuel j on 
Using the solutions of the single-generator dynamic S^ator i at time t under s. For the start-up energy 
programs, one can compute the total generation at each time requirement, we enforce a similar constraint, 
period and under each scenario. We check for optimality 

and, if needed, we repeat the previous process using an ^tfjj -( u h~ u h-0 + su ? 

updated X and p. The violations of the demand and reserve is 
constraints form a subgradient to the dual which can be used 
in the bundle method iterations. 

If a parallel computer is available, the dynamic programs where n ^ the set of facls that may be used to start-up 
of the generating units can be distributed among the pro- generator i. 

cessors. Given that solving the dynamic program of a cycler 20 There might be other constraints that need to be imposed 
requires more calculations than solving the dynamic pro- on fuel allocation. For example, we may want to impose a 
gram of a must-run unit, we distribute the single-generator constraint on the maximum amount of fuel j that can be used 
problems so that the number of cyclers and the number of by a certain generating unit, i, during the planning horizon, 
must-run units are almost the same over all processors. That is, 

25 



Procedure FUEL 



The procedure FUEL is called by the procedure SUCP. 
Procedure FUEL solves a linear program to determine an 
optimal allocation for fuel under certain generation load, 

x, /. The constraint matrix of the linear program is built 30 must be less than the total heat available from fuel j 

iteratively, as shown in FIG. 15. regardless of s. Another constraint might be an upper bound 

The process starts in function block 1501 in which we set 00 lhe amount of fuel i ^ over thc Pining horizon. In 

k«-l, s-2, t-1, i-1, and These variables are the other words ' the value of 

counters used in the construction of the constraint matrix of 35 

the linear program solved to obtain an optimal fuel alloca- V V r» 

tion. In function block 1502, B(k,I(y /f /)) to 1 and B(k,I(y, hh 
<j a(s) )) to -1. Note that I(y M /) is the index of the column in 
the constraint matrix corresponding to the variable y itt f. 

Similarly, l(y ift ° (sy ) is the index of the column in the 4Q must be bounded above by a given constant. An important 
constraint matrix corresponding to y,,,/^. constraint is "fuel j cannot be used on unit i" which can be 
In function block 1503, the value of j is incremented by modeled by setting y />f /-0 for all t and s. To simplify 
1. Then, in decision block 1504, if j=F, the execution is notations, we write fuel constraints (which are assumed to be 
directed to function block 1502; otherwise, the execution is linear) in the form Ay ^b* where A is a constraint matrix of 
directed to function block 1505. In function block 1505, i is 45 appropriate dimension and b A is the corresponding right- 
incremented by 1. Then, in decision block 1506, if i^N, the hand side. 

execution is directed to function block 1502; otherwise, the Note that the fuel usage of a generator must fulfill the 

execution block proceeds to function block 1507. In func- bundle (i.e., non-anticipativity) constraints imposed by the 

tion block 1507, the value of t is incremented by 1. Then, in scenario tree. For example, if Scenario 1 and Scenario 2 

decision block 1508, if tix(s-l), the execution is directed to 50 have the same load between the first period and t, then 

function block 1502; otherwise, the execution proceeds to y f ^-y,. f</ - 2 =0 for all i=l, . . . , n, t=l, . . . x, j«l, . . . , F. The 

function block 1509. In function block 1509, the value of the constraint matrix corresponding to bundle constraints is 

scenario index, s, is increased by 1. Then, in decision block denoted by B. The matrix B is constructed using the vari- 

1510, if s^S, the execution is directed to function block ables a and x returned by TREE and BUNDLE. To do so, 

1502; otherwise, thc execution terminates. 55 I(y iV /) is defined to be the column of B (or A) that 

The cost of generating one MWH of electricity is a corresponds to the variable y itt f. An integer k is used to 

function of the cost of fuel used and its heat content. It is count the number of rows in B. The pseudocode for the 

assumed that we have F different fuels and that the cost of iterative process of building the linear program solved by 

providing a BTU using fuel j at time t is h,/ $/BTU. Note FUEL is as follows (see FIG. 15): 

that s here is a scenario index that reflects load forecasts 60 

provided by the user. When TREE builds on these forecasts • initialization Set k • 1. 
to generate new scenarios, each branch in the resulting tree 

has values of h equal to those of the original scenario that is * General Step. 

used to create this branch. Note that h t f is a function of t p or ^1 ^2 s t 

which allows incorporating the fluctuation in fuel cost, over 65 

the planning horizon, into calculations. Given the generation ^ a11 f=1 T ( f " *)* 

load, x t j* f under scenario s on each unit i at time period t, we 
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-continued 

For all i=l, 

For all y=l F, 

Set B(* t /(^j)) *- 1 and B{k. l(y%*])) «- -1. 
Increment k by 1 . 
Next y 
Next i 
Next ; 
Next s 

At termination, k is the total number of bundle constraints. 
To reduce storage requirement, only the non-zero elements 
of B are stored. 

To utilize fuel optimally, the expected cost of the fuel used 
is minimized by solving the stochastic Linear program 

* N T F F 

*r / /* X Z Z + aeT y + f^zL - ? + * * 

• ^— ' i=i *=i y=i y=i 

E - ^ + £ = «/ -«?/-!) + . 

si/ 1 , 

fly =0, 

y* y. 2 fc0 ' 

where y and y are the vectors of excess and shortage of fuel, 
a and p are penalties on violating the constraint 




and 1 is a vector of ones. The previous stochastic linear 
program is solved using an interior point approach (e.g., I. 
J. Lustig, J. M. Mulvey and T. J. Carpenter, "Formulating 
stochastic programs for interior point methods", Operations 
Research, 39:757-770, 1991). Many commercial packages 
already include convenient routines for such techniques 
(see, for example, the Optimization Subroutine Library of 
IBM). If a parallel computer is available, the process can be 
accelerated as described in I. J. Lustig and E. Romberg, 
"Gigaflops in linear programming", Operations Research 
Letters, 18:157-165, 1996. 

Procedure DP 

The procedure DP is called from the bundle trust region 
method of function block 122 in FIG. 12 which shows the 
procedure SUCP. The flow diagram of the procedure DP is 
shown in FIG. 14. Given the Lagrange multipliers, DP 
determines the status and the load on a generating unit by 
solving a dynamic program. In function block 1401, the 
values of p(s), s=l, . . . , S, are initialized to T. The recursion 
starts at time period T and moves backward until time period 
1 is reached. The time period under investigation for sce- 
nario s is p(s). The scenario index s is set to S. The dynamic 
program of scenario s is solved by moving backward in time 
from period p(s) to t(s). What is being done here is similar 
to the dynamic programming approach used in solving 
stochastic decision processes. See, for example, D. P. 
Bertsekas, Dynamic Programming: Deterministic and Sto- 
chastic Models, Prentice Hall, 1987. 
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In function block 1403, the dynamic program of a(s) is 
solved by moving backward from time period p(s) to time 
period x(s). At this stage both of the dynamic programs of s 
and its predecessor a(s) have been solved. The probability 

5 p(o(s)) of c(s) is incremented by p(s). Then, in function 
block 1404, the value of P(a(s)) is set to x(s)-l. The scenario 
index, s, is decreased by 1 in function block 140S. That is, 
the dynamic program of s is done, and we are ready to move 
on to another scenario s-1. In decision block 1406, a check 
is made to determine if s>=0; e.g., have all scenarios been 
studied? If so, the procedure moves to function block 1407; 
otherwise, the procedure moves to function block 1402 to 
study a new scenario. In function block 1407, the solution of 
the dynamic program has been finished, and we move 
forward to trace the optimal decisions of the dynamic 

15 program. The policy obtained at the end of this stage is the 
optimal policy that is returned to the calling routine. 

Given the data related to a generator, procedure DP solves 
the corresponding single-generator dynamic program and 
returns an optimal schedule, u, for this generator and the 

20 load on it, x, under all scenarios at each period in the 
planning horizon. The stages in our dynamic programming 
formulation are the time periods. Each stage has L y +1 ( . states. 
The first L £ states are those corresponding to the generator 
being on and the second 1 ( . states correspond to the generator 

25 being off. The number of states can be reduced to two in 
some cases as suggested by S. Takriti, J. R. Birge and E. 
Long, ibid. 

Given a scenario tree that is described using x(s) and o(s), 
s-l, . . . , S, DP solves a single-generator dynamic program 
30 by moving backward. A method for organizing the calcula- 
tions and computer storage are described that a high per- 
formance is achieved (see FIGS. 13A and 13B). 

For the following pseudocode of the procedure DP, it is 
assumed that scenarios, s=l, . . . , S, are sorted in the 
35 ascending order of x(s). 

Input, integer T, integer N, integer 1, integer L, integer v, 

integer w, float d(l . . . S,l. . . T)> float p(l . . . S) 
Output, integer u(l . . . S,l . . . T), float x(l . . , S, 1 . . . T) 
Local Variables, integer j, integer s, integer t, integer, P(l . . . 
40 S), integer y(l . . . S) 

Initialization. Allocate enough memory for 

s 

l+£(T-r(5)+l) 

45 

stages. Each stage contains L+l states and each state con- 
stitutes of a real number (cost-to-go) and an integer number 
(the predecessor index). Set 

50 

s 

y*-l + £(r-r(5)+l). 



55 Set p(s)*-T, s=l, . . . , S. For each scenario s»l, . . . , S, set 
y(s)-j and j«-j_(T_ T ( s )+l). 
General Step. 
For each scenario s-S, . . . , 1, 

Starting at t=p(s), move backward using dynamic pro- 
60 gramming until t**t(s). These calculations must be 
performed in the allocated space y(s) to y(s)-(Tt(s)+1). 
Each state must maintain the index of the predecessor 
state that provided an optimal solution for it. 
Starting at t»p(o(s)), move backward using dynamic 
65 programming until t«r(s). These calculations arc per- 
formed in the allocated space y(o(s)) to y(o(s))-(T-T 
(s)+l). Set K(o(s))-r(s)-l. 
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Increment the states of stage y(p(s)) by the values of the 

states in y(s)-(T-t(s)+1) weighted by p^. 
Increment p 0(a) by p 5 . 

Using the initial state of the generator, trace the optimal 
paths to obtain an optimal decision tree. Note that to 5 
find that moving backward requires solving an optimi- 
zation problem at each stage. This problem is that of 
minimizing f//(x,-/)-Vx// subject to g^x./^G^. 
There might be some other bounds on x ( / that can be 
incorporated. For example, if "the load on a generator 10 
cannot exceed (g^,+G ifl )/2 in the first hour after turning 
it on", then we can force this constraint on the dynamic 
program in the row corresponding to "On 1" (see FIG. 
13B). 

To accelerate the calculations, DP is applied to cycling 15 
units only. For must-run units, one can find the optimal loads 
by minimizing the function f(x)-Xx at each time period. In 
general, the peakers can be ignored, u ( -/«0, because they are 
operated in cases of emergency. 

System's Output 20 

Now, the output of the tool according to the invention is 
described. First, the status, u t - /, of a generator, i, is a binary 
variable that indicates whether this generator is on, u t -/«l, 
or off, u^M), at time period t under scenario s. The 2 s 
procedure* guarantees that the resulting values of u,/ are 
feasible; i.e., 

(a) they satisfy the minimum up-time and minimum 
down-time requirements. For example, if the minimum 
up-time, L t -, of a generator, i, is equal to three, then 30 
sequences such as 0110, 0100, and 0101, are infeasible. 
The reason for that is that when the generator is 
switched on; i.e., it has to remain on for three 
periods or more. Note that the sequence 0111 is feasible 

in this example. ^ 

(b) they are consistent between the different scenarios. 
This condition is known as non-anticipativity in sto- 
chastic programming (see, for example, Rockafellar 
and Wets, "Scenarios and policy aggregation in opti- 
mization under uncertainty", Mathematics of Opera- 
tions Research, 16(1): 119-147, 1991). To clarify this 40 
concept, let us pay attention to the beginning of the first 
period, t-1, in our planning horizon. Although there are 
different scenarios that may occur in the future, we 
must make one decision for each generating unit 
regardless of which scenario may occur in the periods 4 5 
to come. This implies that the decision we make must 
be feasible for all scenarios, s=l, . . . , S. To further 
clarify the concept of non-anticipativity, let us discuss 
the set of scenarios presented in FIG. 11. The scenarios 
are organized in a scenario tree. Each branching node 50 
in the scenario tree (which must be at the beginning of 

a period) represents a point in the planning horizon at 
which we can make decisions regarding the status of 
our generators: switch on or switch off a unit. The first 
decision point is the beginning of the first period, t«l. 55 
Note that all scenarios go through this node. Hence, the 
decisions made at this node must be feasible to all 
scenarios (because we do not know which scenario we 
may face in the future). Mathematically speaking, we 
must enforce the constraint u, /ou. Vs and s'=e{l, 2, 
3, 4} for each generator ie{l,' . . . ,'N}. During the first 60 
period, we expect our demand pattern to follow one of 
two paths: either that of Scenario 4 or the path of 
scenarios 1, 2, and 3. Let us assume that our demand is 
moving up during the first period; i.e., we are expecting 
scenario 1, 2, or 3 in the future. When we reach the 65 
beginning of period 2 which is our second decision 
point, we must make decisions that are feasible for 
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scenarios 1, 2, and 3 (because we do not know yet 
which one is going to occur). The mathematical con- 
straint is then u^^u,-/ , Vs and s'e{ 1,2,3} for each 
generator ie{l, . . . ,NT. 
Let us say that during period 2, we observe a demand pattern 
similar to that of scenario 2 (or scenario 3 since they are the 
same so far). At the beginning of period 3, we make a 
decision for each generator such that our decision set is 
feasible for the demand of scenario 2 and that of scenario 3. 
The process continues as long as we have two or more 
scenarios sharing the same demand pattern up until this 
point in time. 

To summarize, the procedure produces an implementable 
set of decisions. That is, decisions that satisfy the math- 
ematical conditions imposed by the structure of the scenario 
tree. The fact that we have a set of scenarios improves the 
robustness of the solution obtained without sacrificing the 
implementability of the decisions among the scenarios. 

Second, if the status of a generator, i, at a certain time 
period, t, and under scenario s is 1, then the procedure 
provides the load (MWH) on this generator during t. The 
load, indicated by x ( . /, satisfies the constraints imposed by 
the operational range of the generator. That is, 
X^/iC,,, Vi e{l, . . . ,N}, Vt e{l, . . . , T}, Vs e{l, . . . 
S}. If u,-/=0, then the load is 0. 

Third, 1 the variables y {>r f of FUEL provides an optimal 
fuel allocation. Given the' heat content, a y - BTU/Ton, of fuel 
j, the expected amount of fuel used by generator i at time t 
is 



Of course, when scenario s is observed, we adopt the fuel 
strategy imposed by the variables y it f. 

Fourth, the procedure also provides the average cost per 
unit of power, "k t s , at each time period under every scenario. 
These values can be used in pricing the power at each time 
period of the planning horizon. Note that as time progresses, 
we obtain a new forecast. The excess capacity in the electric 
system (over the refined forecast) represents the amount of 
power that we should offer to sell in the market. In case of 
having less generation than the new forecasted load, we buy 
the needed power from the open market. 

Fifth, the values of the dual variables of the linear 
program of procedure FUEL and of the stochastic program 
of procedure SUCP provide important sensitivity analysis 
information. For example, the value of \ s indicates the 
change in the generation cost if d/ increases by one unit. The 
dual variables of FUEL denote the economical importance 
of each constraint and how it may affect the value of an 
optimal solution. Note that the dual variables corresponding 
to the bundle constraints in FUEL indicate the importance of 
non-anticipativity for a given scenario and time period. 
Finally, one can use the linear program to study the range in 
which fuel prices can change without changing the optimal 
fuel allocation. 

Note that in reality, the observed scenario is (most likely) 
not a member of the set of the original scenarios used in the 
model. One can adopt the solution of the closest scenario 
and apply it to the real system. 

Commercial packages available in the market deal with 
one forecast and provide a feasible schedule. They do not 
consider the uncertainty in future demands and do not model 
the buying and selling transactions. Also, these packages do 
not incorporate fuel constraints and changes in fuel cost in 
their models. On the other hand, our model considers 
uncertainty by using a set of scenarios and produces prices 
at which power may be traded without violating fuel con- 
straints. 
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While the invention has been described in terms of a 
single preferred embodiment, those skilled in the art will 
recognize that the invention can be practiced with modifi- 
cation within the spirit and scope of the appended claims. 

Having thus described my invention, what I claim as new 
and desire to secure by Letters Patent is as follows: 

1. In a risk management system that schedules generating 
units of an electric utility while taking into consideration 
power trading with other utilities and a stochastic load on the 
utility system, a computer implemented method comprising 
the steps of: 

inputting load forecast, trading forecast and generator and 
fuel properties; 

generating load scenarios as a scenario tree that approxi- 
mates future uncertainties in load and trades; 

computing an initial approximate value for the cost per 
British Thermal Unit (BTU); 

obtaining a load on each generator under any scenario 
and, given the load on each generator, determining an 
optimal fuel allocation to obtain a cost function of 20 
generating a Megawatt Hour (MWH) of power on each 
generator; 

solving a stochastic unit-commitment model to obtain a 
generation of each unit at each time period under 
different scenarios which can be used to compute 
heat-input needed for each unit throughout a planning 
horizon; 

solving a linear program in which decision variables are 
amount of fuel burnt on each unit at each time period, 
wherein decision variables have a component for each 
possible fuel in the model and the linear program is 
solved subject to a constraint that fuel fed to a generator 
must be sufficient to generate the heat needed; 

determining if a current fuel allocation is the same as 
before and, if not, computing a new cost function for 
each generator for all fuels fed to a given generator, the 
new cost function being used in resolving unit com- 
mitment; 

repeating the steps of solving the stochastic unit- 
commitment model until the fuel mixtures and the load 
on the generators approach an optimal solution; and 

outputting fuel consumption, generation requirement and 
load for each generator of the utility at each period and 
the fuels to be used. 

2. The computer implemented method of claim 1 wherein 
a generator may use more than one fuel at any time period, 
the step of computing an initial approximate value for the 
cost per BTU being performed by using the fuel with the 
least cost. 

3. The computer implemented method of claim 1 wherein 
the step of obtaining a load on each generator is performed 50 
by the steps of: 

constructing cost functions using an average fuel cost; 

for each scenario and period in a planning horizon, 
ranking the generators in an ascending order of average 
marginal cost at a maximum production level, selecting 55 
enough generators so that the load and spinning-reserve 
constraints are met, and solving an economic power 
dispatch; 

determining a feasible solution for a stochastic unit com- 
mitment model; and 
solving an economic power dispatch. 

4. The computer implemented method of claim 1 wherein 
a generator may use more than one fuel at any time period, 
the step of computing an initial approximate value for the 
cost per BTU comprises the steps of: 

approximating the cost per BTU by using a fuel with the 
least cost; 



solving a unit commitment problem to determine the 

amount of power generated by each generating unit; 
finding an optimal fuel mixture for each generator at each 
time period and under each scenario using linear pro- 
gramming; and 
repeating the previous two steps as long as the fuel 
mixture obtained from the linear programming solution 
changes between iterations. 

5. The computer implemented method of claim 1 wherein 
a user of the risk management system may provide a 
scenario tree which is processed to define bundle 
constraints, two scenarios being members of the same 
bundle at time t if their loads are the same up until time t. 

6. The computer implemented method of claim 5 wherein 
the processing of the scenario tree to define bundle con- 
straints comprises the steps of: 

setting a predecessor of Scenario 1 to 1, a split time of 

Scenario 1 to 1, and a scenario index, s, to 2; 
comparing a load of Scenario s to a load of Scenario 1 to 
determine a split time of s as a first point in time at 
which the load of Scenario s is different from that of 
Scenario 1 and, at this point, setting a predecessor of s 
to 1; 

introducing a new variable, s', and initializing s' to 2; 
comparing the load of s to all other scenarios s'-l, . . . , 
s-1, in order to determine a split time and predecessor 
of s; 

determining whether s' is less than s and, if so, incre- 
menting s by 1, but otherwise, comparing split times of 
s and s'; 

if split times of s and s* are not equal, indicating that s' 
should not be used as a predecessor to s, incrementing 
s f by 1, but otherwise, comparing the loads of s and s' 
to obtain a new split time for s and setting the prede- 
cessor of Scenario s to s 1 before incrementing s* by 1; 
when all the scenarios in the scenario tree are scanned, 
checking s against a total number of scenarios, S, and 
if s=S, terminating the execution but otherwise, the 
execution is continued until all scenarios, S, have been 
processed. 

7. In a risk management system that schedules generating 
units of an electric utility while taking into consideration 
power trading with other utilities and a load on the utility 
system, a computer implemented method comprising the 

45 steps of: 

inputting factor data comprising load forecast, trading 
forecast, spot market estimates, and generator proper- 
ties and fuel properties, costs and constraints; 
generating load scenarios that approximate future uncer- 
tainties in at least some of the factor data; 
obtaining a load on each generator under at least one of 
said scenarios and determining an optimal fuel alloca- 
tion to obtain a cost function of generating power on 
each generator; 

solving a plurality of commitment models to obtain a 
generation of each unit at each time period under 
different scenarios which can be used to compute 
energy input needed for each generator unit throughout 
a planning horizon; 
selecting an optimal fuel mixture and a load on the 

generators based on the models; and 
outputting fuel consumption, generation requirement and 
load for each generator of the utility at each period and 
the fuels to be used and outputting average production 
65 costs at each period under each scenario. 
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